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1. Introduction 

One of the hardest problems in Lattice QCD is the calculation of disconnected quark loops 
and all-to-all quark propagators [|l|]. Statistical noise methods can, in principle, make any matrix 
element available, but the computer expense can be prohibitive. This is especially true at low quark 
mass, where error bars are greatest. Noise subtraction methods, such as the perturbative subtrac- 
tion method |^, which improve the signal by a suppression of off-diagonal noise, are essential 
for efficiency. However, perturbative subtraction is ineffective at low quark mass. We propose a 
subtraction method, which we call "eigenspectrum subtraction", which uses low eigenmode infor- 
mation to suppress the statistical noise at low quark mass. Such eigenmode information is readily 
available from fermion deflation algorithms, such as GMRES-DR for non-hermitian systems or 
LAN-DR [Q] for hermitian ones. 

2. Method 

2.1 Introduction 

Our noise methods will be utilizing real Z(2) noise. These are vectors made up of statistically 
random I's and — I's. A useful property of these Z(2) noise vectors is 

The quark propagator problem can be formed by solving for the solution vector of a linear 
system made with the necessary action, M, and any given source vector, b. 

Mx = b 

, (2.2) 

For our noise algorithms we will be using many Z(2) noises as our source vectors so we get 

^=lMr^^^ (2.3) 
This formulation is done with a Wilson action of the form 

M=\-kD (2.4) 

For the disconnected quark loops and all-to-all quark propagators we will, as a test, form 
Tr(M^'), although the method can be used for any quark matrix element. 

2.2 Non-Subtraction 

First let us consider the the basic setup. By utilizing the aforementioned properties we can 
compute the trace of the inverse of a matrix. 
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This method allows for the trace computation to be done by solving many linear systems. 
The error from this method is introduced by the off-diagonal elements, that is; the less diagonally 
dominant the matrix is, the greater the error that is introduced. In order to minimize this effect we 
must form a matrix that approximates the off-diagonal elements of M^^. We will call such a matrix 
M^^. We will use this matrix to subtract out the off-diagonal elements. 

2.3 Perturbative Subtraction 

A useful method is to use the perturbative approach to form the approximation matrix[^. That 
is to form 

Mpelt = 1 + kD+ {xDf + {xDf + {xDf (2.6) 

Since M is in a form that is conducive to small parameter expansions we can subtract off the 
necessary, approximate, off-diagonal elements. Since Mp}rt is only an approximation it does have 
on-diagonal elements and thus the trace of Mp}^, must be added back. 

Tr(M-i) = jim zJ)£(M,^^ _M^;jj4«) +Tr (M,;,) (2.7) 



This method produces significant improvements to the "Non-Subtraction" method, but with a 
h. Since Mp}^^ is a small parame 
relate to "large" quark mass values. 



catch. Since Mp}^^ is a small parameter expansion, it is most effective for small values of k, which 



2.4 Eigenspectrum Subtraction 

In order to allow for small values of the hopping parameter we must formulate a new ap- 
proximation matrix. Our new method allows for an eigenspectrum formulation by utilizing the 
eigenvectors, both left and right, from the non-hermitian Wilson matrix, M. 

Right eigenvectors are formed as 

m4'''=A('?)4'^ (2.8) 
where left eigenvectors are formed like so 
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4«)'m = 4'^'a(«) (2.9) 

It is useful to note that the eigenvalues of the left and right eigen-system problem are the same. 
That is to say for every right eigenvector, there exists a left eigenvector. Left and right eigenvectors 
are orthogonal. (For the Wilson case, the two are simply related by 75.) Thus any matrix, mxm 
can be formed via an eienspectrum formulation by 

M= f;4''^^^^4^^' (2.10) 

By using the orthogonality of this form and only using the Q smallest eigenvalues we form our 
approximation. 



^..-L^(,)-«-L ^^^^^ 
where e]^^ ■e'f^ = 

The smallest eigenvalues of M are used as we expect their contributions will be the greatest. By 
utihzing this eigenspectrum formulation, we are better able to represent the off-diagonal elements 
for small quark masses. Since at small quark masses the eigenvalue spectrum shows many small 
eigenvalues, we expect that we will get increasingly better results in a the range of Kcru- 

As before we can form the trace as 

UnUke the Perturbative method, however, the Tr ( M^J ) can be formed with extreme ease. 




(2.13) 



The ease of this method does not end here. Since no part of Mgig is ever actually formed, 
but rather the right/left eigenvectors are, the number of matrix/vector multiplications is drastically 
minimized when compared to the Perturbative Expansion Method. 



3. Perturbative Subtraction and Eigenspectrum Subtraction 

In an attempt to combine both method in hopes to yield a "hybrid" method we developed 
something new. To naively apply both methods at once would form an approximation matrix of the 
form 

M~^- =M~^ +M~^ (3 I) 

^"naive ^"pert ^ ^"eig ■ \~>-'^J 

The issue with this naive formulation is the the approximate off-diagonal information would 
be removed twice over, introducing error. In order to correct for this issue we need to remove the 
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low eigenspectrum information from MpJrt and replace it with the low eigenspectrum information 
ofAf-i. This new formulation would provide a perturbative approximation, but where the low 
eigenmode information of the approximate perturbative matrix is replaced with the low eigenmode 
information of M. This should allow us to approximate for both large and small values of K. To 
subtract the low eigenmode information from Mpen we form it as 



^ 1 (n\ (n\T 

1 ^ (3.2) 



Where ^ = .?)m;;,4^) 

T] is an eigenvalue-like term that approximates the eigenvalues of M~g„. With this formulation 
our new trace comes to 



n i 

(3.3) 



= li^i ^ 1 4f ^ - M;/,_, . - )z W + 1 ^ + Tr {Mp},) - £ -L 

n i q q 'I 

4. Larger Problems 

The Wilson action can be formed as 

"--(it) 

allowing for the linear system to be solved as 

Due to the form of the matrix, a "reduced matrix" can be formed allowing for the problem to 
be cut in half. 

^reduced ~ o ^eo^^oe (4.3) 

The reduced linear system then becomes 

MreducedXe = A^e + r^eoK (4.4) 

and Xo can then be directly computed by 

Xo = bo + KHoeXe (4.5) 
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When programming for Lattice QCD, the reduced problem is more efficient. The Eigenspec- 
trum Subtraction information can be extracted from the reduced system in the Wilson case. The 
eigenspectrum of the reduced matrix will be formed in the usual way. 

Mreducedxi'^ =X^'l^xi''^ (4.6) 

The eigenvalue of the reduced matrix (X), can be related to the eigenvalue of the full matrix 

(A). 

A = 1 ± \Jl-K^X (4.7) 

The eigenvectors, likewise, can be be related. Xe, the eigenvector of the reduced matrix, is also 
the even part of the eigenvector of the full matrix. 



Since the reduced matrix can be used to form the eigenspectrum information of the full matrix, 
which is readily available by linear solver algorithms with deflation, forming the full matrix is never 
needed. This is untrue for the Perturbative Subtraction method. The low eigenspectrum information 
is related to the "negative sign" of X and Xg. 



5. Tests and Results 

As a test, we applied the method to an 8^ Wilson matrix, M, using a parallel version of MAT- 
LAB. The Kcrit = 0.15701 value was determined. 

Since we are only interested in comparing size of the error bars from different methods, and 
thus each methods effectiveness various values of K, we zero out all the values themselves in the 
figures and only show the error bars. 

The trials consisted of runs where 100 real Z(2) noises were used for 20 different times, each 
time with a different random seed, to define error bars. We show the results after 1, 10, 20,..., 
100 noises, or iterations. In order that error bars do not overlap, the results are separated from 
one another by a small value on the noise axis. (The non-subtracted results are located exactly at 
positions 1,10, 20, 30,..., 100). 

The notation is as follows. 



NS 


Non-Subtracted 


PE 


^th Order Perturbative Subtraction 


eev 


Eigenspectrum Subtraction with Q eigenvectors 


PEc+2 ev 


Perturbative Subtraction corrected with Q eigenvectors and 
Eigenspectrum Subtraction with Q eigenvectors 
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(b) k: = 0.1560 
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Figure 2: Plot of the low eigenspectrum of the 8^ Wilson lattice at Kcr = 0.15701. 



6. Conclusions 

At a value of K near K^rit our method shows evidence of improvement over the 4th order 
Perturbative Subtraction Method. By combining the two methods, and removing the low eigen- 
information from the Perturbative Subtraction Method, we see an additional increased improve- 
ment, albeit small. As K deviates from Kent we see results on par with Perturbative Subtraction. 

We would expect our method to become more efficient for larger matrices with smaller eigen- 
values, but we need more study. The combined effects of Perturbative Subtraction and Eigen- 
spectrum Subtraction are still being studied as well. Development in FORTRAN is in progress to 
facilitate the larger matrices. 
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